Isotopic chains around oxygen from evolved chiral two- and three-nucleon interactions 
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We extend the formalism of self-consistent Green's function theory to include three-body interactions and 
apply it to isotopic chains around oxygen for the first time. The third-order algebraic diagrammatic construction 
[ADC(3)] equations for two-body Hamiltonians can be exploited upon defining system-dependent one- and 
two-body interactions coming from the three-body force, and correspondingly dropping interaction reducible 
diagrams. This goes beyond the standard normal ordering truncations recently used in ab-inito studies. The 
Koltun sum rule for the total binding energy acquires a correction due to the added three-body interaction. 
This formalism is then applied to study chiral two-nucleon (2N) and three-nucleon forces (3NF) evolved to 
low momentum cutoffs. We find that these interactions reproduce the binding energies of nitrogen, oxygen and 
fluorine isotopes to great accuracy, providing clear indication of the predictive power of this approach. All three 
neutron driplines are correctly predicted when full 3NF are included. The formalism introduced also allows to 
calculate form factors for nucleon transfer on doubly magic systems. 

PACS numbers: 21.10.-k, 21.30.Fe, 21.60.De 
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Introduction. - The ultimate goal of ab-initio nuclear the- 
ory is to achieve accurate predictions of nuclear properties that 
are consistent with the underlying theory of QCD and are pa- 
rameter free as much as possible. Advancing on this problem 
is presently of primary importance in the mid mass region of 
the nuclear chart, in response to significant advances in the 
discovery of new nuclides at radioactive isotope facilities [1]. 
Moreover, parameter free predictions would help reducing un- 
certainties in our knowledge of dripline isotopes that are cur- 
rently beyond experimental reach [2]. The requirement of ac- 
curacy poses strong requirements on both the knowledge of 
the nuclear interaction and the many-body techniques used in 
the calculations. 

Breakthroughs over the last decade were possible due to 
the introduction of many-body methods that scale gently with 
increasing particle number and therefore can reach systems 
well above the p shell. Self consistent Green's function the- 
ory (SCGF) [3, 4], coupled cluster (CC) [5-7] and in-medium 
similarity renormalization group (IM-SRG) [8, 9] have been 
employed in ab-initio calculations of doubly closed shell nu- 
clei with masses up to A~60. For open-shells, semi-magic 
isotopes can be calculated by breaking particle conservation 
symmetry and reformulating theories in terms of Hartree-Fock 
Bogolioubov reference states as done in Grokov theory [10] 
and in IM-SRG. Calculations based on IM-SRG have been 
performed for ground state energies [9]. On the other hand, 
the state-of-the-art SCGF theory not only can be extended to 
open shells [11] but it gives access to a wealth of nuclear struc- 
ture information. This includes the addition or removal of one 
or two nucleons to/from the calculated ground states [12, 13] 
and direct link to microscopic optical potentials [14]. 

In this Letter, we extend the scope of SCGF to include 
three-nucleon forces (3NFs) in finite nuclei. We define density 
dependent one- and two-body forces derived from the 3NF 
part of the hamiltonian and work out the first order correction 
to the Koltun sum rule to obtain binding energies. The method 



is applied to the isotopic chains of oxygen, nitrogen and fluo- 
rine, as well as spectra of single-neutron states in the sd shell. 
This opens the possibility of probing modern realistic nuclear 
interactions on a wide range of experimental data, including 
excitation spectra, the evolution of shell closures, and the po- 
sition of driplines. 

Formalism. We employ Green's function (or propagator) 
theory, where the object of interest is the single particle prop- 
agator [15]. 
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where greek indices afi,..., label a complete orthonormal ba- 
sis set and s^ +l = (£^ +1 - E%) and e^ 1 = (EJ* - E^ 1 ) 
are the nucleon addition and separation energies, respectively. 
In Eq. (1), |^ +1 >, |^ _1 > are the eigenstates and E* + \ E^ 1 
are the eigenenergies of the (A + l)-nucleon system. Hence, 
gapiw) describes the exact propagation of a single-nucleon or 
hole excitation through the system. From Eq. (1) we also ex- 
tract the one-body reduced density matrix, 
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where the integration contour CT is taken on the upper half of 
the imaginary plane. 

We start our calculations with the intrinsic Hamiltonian 
H(A) = H - T CM ,(A) = (7(A) + V(A) + W in which the kinetic 
energy of the center of mass (c.o.m.) has been subtracted and 
we put in evidence the dependence on the number of nucle- 
ons A. The terms U, V and W collect all the the one-, two- 
and three-nucleon contributions, respectively. Based on this, 
we define system dependent one- and two-body effective in- 
teractions obtained by contraction with the correlated density 
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matrix, Eq. (1), 
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where all matrix elements are properly antisymmetrized and 
summation over repeated indices are implied here and in the 
following. The resulting Hamiltonian, H — U + V+W, can be 
proved to lead to the same Green function (1) for the original 
hamiltonian with the caveat that only interaction irreducible 
must be retained in their diagrammatic expansion [ 1 6] 1 . Equa- 
tions (3) and (4) generalize the idea of normal ordering of the 
Hamiltonian to fully correlated densities. In the following we 
retain only the U and V terms and discard diagrams with ex- 
plicit interaction irreducible 3NF terms. This approximation 
has been proven to be negligible for the present hamiltoni- 
ans in [17, 18] 1 . The single-particle propagator g a p{<x>) can 
then be calculated for finite closed-shell nuclei by exploiting 
the effective one- and two- body interactions with the already 
available two-body formalism. 

We first solve the spherical Hartree-Fock (HF) equations 
for the full Hamiltonian within the given model space. The re- 
sulting propagator, g^(oj), is then used as a reference state to 
calculate the energy-dependent contribution to the self-energy. 
We employ the third-order algebraic diagrammatic construe - 
tion [ADC(3)] method [19, 20] and write the self-energy as 
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where M (N) are interaction matrices in the 2hlp (2plK) space 
and C (D) are the corresponding coupling strength to the sin- 
gle particle states. In the ADC(3), these matrices are con- 
structed to guarantee that all diagrams up to third order are 
included in Eq.(5). We note that the ADC(«) approach de- 
fines a hierarchy of truncation schemes for increasing order 
« to guide systematic improvements of the method. The cor- 
related propagator g a p{oS) is finally obtained by solving the 
Dyson equation, 
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which is diagonalized using a Lancsoz algorithm as explained 
extensively in [21, 22]. Note that we employ the scO approxi- 
mation of Refs. [11, 22] in which only the EJaf) contribution 
of Eq. (5) depends on the reference states g^{a>). This im- 
plies the iterative solutions of Eq. (6) to evaluate ~L^=U a p in 
terms of the final correlated density matrix, Eq. (3). 



1 A diagram is said to be interaction reducible if it can be factorized in two 
lower-order diagrams by cutting an interaction vertex or, equvalently, if 
there exists a group of lines (interacting or not) that leave an interaction 
vertex and eventually all return to it. 

2 This would be the analogous to the N02B approximation of Ref. [18] but, 
here, not tight to the choice of any reference state. 



In presence of 3NFs, the ground state energy can still be 
inferred from the Koltun sum rule (SR) that now acquires a 
correction: 

(7) 

Eq. (7) is still an exact equation. However, it requires to eval- 
uate the expectation value of the 3NF part of the hamiltonian 
(VPqIW'II'q), with an accuracy comparable to the many -body 
approximation in use. We calculate this correction at first or- 
der in W using fully correlated propagators, 
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that implicitly includes relevant higher order terms from stan- 
dard many-body perturbation theory. We found that it is 
mandatory to use the fully dressed propagators, i.e. solutions 
of the Eq. (6), but this is sufficient to account for all relevant 
contributions. The next order correction is given by 
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where AT is the two-body density matrix after subtraction of 
the zeroth-order contribution coming from two fully corre- 
lated but non-interacting nucleons, to avoid double counting 
Eq. (8). We estimated this using in Tamn-Dancoff approxima- 
tion (TDA) [23] and found its contribution to be small com- 
pared to our estimated errors, as discussed below. 

The binding energy and spectra of neighboring even-odd 
isotopes are extracted from the poles of propagator (1), how- 
ever this requires a proper correction to account for the vari- 
ation in the kinetic energy of the c.o.m. motion with chang- 
ing A. To extract the energy of a system with mass A + 1, we 
recalculate g a p{co) for the doubly closed-subshell A-nucleon 
system but with a H(A +1) corrected hamiltonian. We then 
obtain: 

E A±1 = +e A±1 [H(A + 1)] + E A [H(A + 1)] , (10) 

were we made explicit the dependence on c.o.m. correction 
for the Hamiltonian and E A [H(A + 1)] is obtained from the 
corrected Koltun SR. 

Results. We perform calculations using chiral effective field 
theory (EFT) two-nucleon (2N) and 3NFs evolved to low mo- 
mentum scales by using using free-space similarity renormal- 
ization group (SRG) [24, 25]. The original 2N interaction is 
N 3 LO with cutoff A 2/v =500 MeV [26, 27] while for the 3NF 
we chose a N 2 LO in a local form [28] with a reduced cut- 
off of A3at=400 MeV and low-energy constants cq=-0.2 and 
ce =0.098 refitted to reproduce the 4 He binding energies, as 
discussed in Ref. [29]. The SRG evolution on the sole 2N- 
N 3 LO interaction already generates 3N interactions, which 
we will refer hereafter to as "induced" 3NFs. When the pre- 
existing 3N-N 2 LO interaction is included in the SRG evoled 
interaction will be referred to as "full" 3NFs. The two-pion 
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FIG. 1. (Color online) Convergence of 16 and 24 as a func- 
tion of increasing size of the model space for Tiuj=2A MeV and 
'lsffG-2-0 frrr 1 . Left: Binding energies from the corrected Koltun 
sum rule, Eq. (7), as obtained from the induced and full interac- 
tions. Right. Expectation value (Y^WI 1 ?^ } obtained at first order 
only, Eq. (8), (full lines) and with correction from two-body TDA 
ladders, Eq. (9) (dashed lines). 



exchange diagrams that incorporate physics from the Fujita- 
Miyazawa 3NF appear in the chiral 3N-N 2 LO force and there- 
fore their effects are incorporated only in the full Hamiltonian. 
Calculations were performed in model spaces up to 12 har- 
monic oscillator shells [N max = max (2« + /) = 11], includ- 
ing all 2N matrix elements and limiting 3NF ones to config- 
urations with N1+N2+N3 < NzNF,max-^- Except when indi- 
cated, a SRG cutoff /i=2.0 firr 1 and harmonic oscillator fre- 
quency hco-24 MeV were used. 

Figure 1 shows the convergence pattern of the total bind- 
ing energy of le O and 24 as a function of the model space 
size. The convergence is optimized by the choice of the cho- 
sen oscillator frequency which is close to the minimum of the 
fiD dependence for the present interaction [9]. The stagger- 
ing between adjacent results is due to the particular trunca- 
tion of 3NF matrix elements and the alternate parities of har- 
monic oscillator shells. Separate exponential fits to the cal- 
culated 24 energies, for N max either even or odd, differ by 
100 keV and are within 600 keV of the N max -l 1 result, with 
an overall error of 0.5% of the total biniding energy. The right 
panel of Fig. 1 demonstrates the similar convergence of the 
(T^IWIT^) expectation values. The contribution of (W TDA ), 
Eq. (9), is never bigger than 300 keV and represent a small 
correction to the Koltun SR. A proper study of the contribu- 
tions to < W > will be addressed in a forthcoming publication. 
For the present purposes, we sum the above uncertainties to 
make a conservative estimate for our convergence error of 1% 
for the calculated total binding energies. 

The differences between calculated binging energies the ex- 
periment, are demonstrated in Fig. 2 for three different com- 
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FIG. 2. (Color online) Differences between calculated and experi- 
mental ground state energies of oxygen isotopes for the full interac- 
tion with different values of Tia> and Asrg- 



binations of the tia) and Asrg- Changes with the variation of 
these parameters are rather small and increase only slightly 
going to larger masses. For 24 we find about 450 keV vari- 
ations for frco between 20 and 24 MeV. This is 0.3% of the 
total energy, which is well within our estimated error from 
convergence. Changes from Asrg m the range 1.88-2.0 firr 1 
are also small and in agreement with the 1% estimated from 
larger ranges of the cutoff [9]. 

From previous studies based on the ADC(3) method, we 
expect accuracy of the many-body truncation scheme of about 
1% [30, 31]. The extrapolated ie O ground state (Fig. 1) is 
over bound at -130.8(1) MeV but in close agreement with 
the -130.5(1) MeV obtained from IM-SRG [9], giving further 
confirmation of the accuracy achieved by different many -body 
methods. Binding energies of 15 and 23 can be obtained in 
two different ways, from either neutron addition or removal 
on neighboring sub-shell closures. Results differ by at most 
400 keV, again within the estimated uncertainty of our many- 
body truncation scheme. The correction (10) is crucial to ob- 
tain this agreement. For hu>=24- MeV and Asrg-2.0 fin -1 , 
the discrepancy in 15 ( 23 0) is 1.65 MeV (1.05MeV) when 
neglecting the changes in kinetic energy of the c.o.m. but it 
reduces to only 190 keV (20 keV) when this is accounted for. 
This gives us confidence on the procedure to separate the ef- 
fects of the center of mass motion and an additional confirma- 
tion of the quality of the method. 

Fig. 2 also gives a first remarkable demonstration of the 
predictive power of chiral 2N+3N interactions: accounting 
for the precision of our many-body approach and dependence 
over larger ranges of cutoffs we expect an accuracy of «3% 
on binding energies and all calculated values agree with the 
experiment within this limits. Note that the interaction em- 
ployed were only constrained by 2N and 3 H/ 4 He data. 

Figure 3 collects our results for the oxygen isotopic chain, 
as calculated for tia>-24 MeV and Asrg=2-0 fin -1 . The top 
panel demonstrates the evolution of neutron single particle 
spectrum in the sd shell. Induced 3NFs reproduce the overall 
trend but predict a bound 1J3/2 orbit for the neutron rich iso- 
topes. The main consequence of adding pre-existing 3NFs — 
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FIG. 3. (Color online) Top. Evolution of single-particle energies for 
neutron addition and removal around sub-shell closures of oxygen 
isotopes. Bottom. Binding energies obtained from the Koltun SR and 
the poles of propagator (1), compared to experiment (bars) [32, 33]. 
All points are corrected for the kinetic energy of the c.o.m. motion. 
For all lines, red squares (blue dots) refer to induced (full) 3NFs. 



magic odd-even isotopes of nitrogen and fluorine. Induced 
3NF forces consistently under bind these isotopes and even 
predict a 27 N close in energy to 23 N. This is fully corrected by 
full 3NFs that strongly binds 23 N with respect to 27 N, in accor- 
dance with the experimentally observed dripline. The repul- 
sive effects of filling the 1J3/2 is also observed in 29 F. However, 
the inclusion of an extra proton provides enough extra binding 
to keep the latter isotope bound by about 700 keV with respect 
to 25 F, in much better agreement with the experimental value 
of 1.47 MeV. The induced interaction alone would overesti- 
mate this binding and pre-existing 3NFs are fundamental in 
achieving the correct balancing between the attraction gener- 
ated by the extra proton and the repulsion due to the filling of 
the neutron sd shell. 
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FIG. 4. (Color online) Binding energies of odd-even nitrogen and 
fluorine isotopes calculated for induced (red squares) and full (green 
dots) interactions. Experimental data are from [33]. 



the full Hamiltonian — is to raise this last orbit above the con- 
tinuum threshold and confirms the increasing repulsive ef- 
fects of the two-pion exchange Fujita-Miyazawa interaction 
on this orbits, as the neutron sd shell is filled [34]. Instead, 
the ^5/2 quasiparticle states are lowered by about 1 MeV on 
average, providing extra binding through the Koltun SR for- 
mula (7). The consequences of this trends are demonstrated 
by the calculated ground state energies shown in the bottom 
panel: the induced hamiltonian systematically under binds the 
whole isotopic chain, and confirms earlier predictions based 
solely on the original 2N-N3LO interaction [35]. The dripline 
is also erroneously placed at 28 because of the lack of re- 
pulsion in the ^3/2 orbit. On the other hand, contributions 
from pre-existing 3NFs are substantial and increase with the 
mass number up to 24 0, when the unbound d3/2 orbit starts be- 
ing filled. As a result, the full Hamiltonian nicely reproduces 
both the experimental ground state energies and the observed 
dripline at 24 [36]. Our result suggest a ground state reso- 
nance for 2s O unbound by 5.2 MeV with respect to 24 0. How- 
ever this estimate is likely to be affected the presence of the 
continuum which is important for this nucleus but neglected 
in the present work. 

The same effects are demonstrated in Fig. 4 for the semi- 



In conclusion, we have considered the extension of the 
SCGF method to include three-body hamiltonians. By prop- 
erly defining system dependent effective one- and two-body 
interactions that include the relevant contribution form 3NFs, 
calculations can be performed with formalisms already ex- 
isting for two-body Hamiltonians. This approach, however, 
goes beyond usual truncations based on normal ordering of 
the Hamiltonian and employs fully correlated densities instead 
of unperturbed reference states. We applied this approach for 
the first time to study SRG-evolved chiral 2N and 3N inter- 
actions on the isotopic chains of nitrogen, oxygen and fluo- 
rine. We find that chiral 3NF at N2LO are crucial in predicting 
the binding energies of these isotopes and they reproduce the 
correct behaviour at the neutron driplines for all three cases. 
Within the estimated errors due to the many-body techniques 
and the dependence on the SRG evolutions, we find a remark- 
able agreement between our calculations and the experimental 
energies along all three isotopic chains. 

Recent results [11] clearly show that state of the art SCGF 
methods can be straightforwardly extended to the correspond- 
ing Gorkov formalism for open shells, which is now under- 
way. This would not only allows direct calculations of semi- 
magic even-even isotopes with analogous quality as above but 
would also allow extracting a wealth of information on neigh- 
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bor isotopes reachable by transfer of one or two nucleons. 
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